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ABSTRACT 

We have analyzed a large sample of numerically simulated clusters to demonstrate the adverse 
effects resulting from use of X-ray fitted /3-model parameters with Sunyaev-Zeldovich effect (SZE) 
data. There is a fundamental incompatibility between /3 fits to X-ray surface brightness profiles and 
those done with SZE profiles. Since observational SZE radial profiles are in short supply, the X- 
ray parameters are often used in SZE analysis. We show that this leads to biased estimates of the 
integrated Compton y-parameter inside rsoo calculated from clusters. We suggest a simple correction 
of the method, using a non-isothermal /3-model modified by a universal temperature profile, which 
brings these calculated quantities into closer agreement with the true values. 

Subject headings: galaxies:clusters:general-cosmology:observations-hydrodynamics- 

methods:numerical-cosmology:cosmic microwave background 



1. INTRODUCTION 

The hot gas in clusters of galaxies is responsible for in- 
verse Compton scattering cosmic microwave background 
(CMB) photons as they travel through the intraclus- 
ter medium (ICM). This results in a spectral distor- 
tion of the CMB at the location of clusters o n the 
sky, referred to as the lSunvaev fc Zeldovic h (1972) effect 
(SZE). This distortion is characterized by a low frequency 
(<218 GHz) decrement, and higher frequency (>218 
GHz) increment in the CMB intensity (jCarlstrom et alJ 
2002). TheX -ray emission in clusters consists of thermal 
bremsstrahlung and line emission from the same highly 
ionized plasma that scatters the CMB. 

High resolution X-ray or SZE observations of clusters 
coupled with assumptions about the gas distribution lead 
to estimates of the gas mass in the cluster dark mat- 
ter potential well. The elect ron number density is often 
assum ed to fit a /3 model (jCavaliere fc Fusco-Femia^nol 

n e (r)=n e0 11+ (-) J . (1) 

In the above relation, n eu is the central density normal- 
ization and r c indicates the fitted parameter referred to 
as the core radius. Fitting an observed X-ray or SZE 
profile to these projected /3 model X-ray surface bright- 
ness and SZE y parameter distributions results in a de- 
scription of the density distribution, which can be inte- 
grated to obtain the gas mass. The difference in depen- 
dence on gas density and temperature of X-ray emissiv- 
ity and the SZE y parameter makes the combination of 
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these two methods of observation potentially very power- 
ful. Because of this difference, the observability of clus- 
ters via each method is affected differently by the im- 
pact of physics in cluster cores including radiative cool- 
ing and feedback mechanisms, as well as the transient 
boosting of surface brightness and spectral temperature 
generated during m erging events (|Roettiger et al.lll996t 
iMotl fc Burnsll2005[ ) . These two methods not only select 
a different sample of clusters, but combined SZE/X-ray 
observations of individual clusters allow one to extract 
the density and temperature of the gas without relying 
on X-ray spectral temperatures. 

1.1. Isothermal Beta Models 

Under the assumption that the gas in clusters is 
isothermal, one can fit an isothermal /3-model to the data 
in order to deduce the density profile. To generate the 
projected X-ray surface brightness profile, we integrate 
the /3-model density distribution 

Sx = ^ +z)i J n e (r)n H (r)A(T)dl, (2) 

where in the bremsstrahlung limit, A(T) oc T 1 / 2 , and in 
a fixed X-ray band is mo re weakly dependent on temper- 
ature (|Mohr et al.| [l999). This integration results in 

Sx(b) = S X0 I 1+ f-J J , (3) 

where Sxo is the fitted central X-ray surface brightness of 
the model, and b indicates the projected radius. Similarly 
for the SZE, a /3-model density profile can be integrated 

/kbT 
o- T n e {r) g dZ, (4) 
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which results in a projected radial distribution of the 
Compton y parameter 

(I_M 
1+ (£) ) • (5) 

For very hot clusters (T > 10 keV), relativistic 
corrections must be included to the SZE integral 
(|Itoh k. Nozawal I2004D . While these fits result in a 
description of the cluster density profile, we must be 
aware that that description is only approximate, since 
it has been shown both in simulations a nd observations 
(|Loken et al.ll2002tlVikhlinin et al.ll2005f) that many clus- 
ters show a radial dependence of temperature. That 
means that when using isothermal models (the above 
equations) to fit the data, we should expect error to be 
introduced in the derived quantities. 

1.2. SZE/X-ray Derived Quantities 

Recent studies have used the values of the /3-model pa- 
rameters determined from the X-ray surface brightness 
profiles combined with SZE cluster observati ons to deter- 
mine the value of the Hubble constant (Ho) (|Reese et al.1 
l2002tlBonamente et al.ll2006l) and the c luster gas fraction 
(|Jov et alJ 120011: lLaRoaue et al.ll2006D . While joint fits 
of X-ray and SZE interferometric data are used to deter- 
mine the /3-model parameters in these studies, the X-ray 
data drives the fit, since the SZE data currently lacks 
the resolution (and interferometri c U-V plane coverag e) 
to constrain the parameters well (|LaRoque et al.1 12006). 

An additional calculation can be performed using com- 
bined SZE/X-ray data from clusters. It is expected both 
from analytic arguments and numerical simulations that 
the integrated SZE signal, as a measure of the total ICM 
pressure in cluste rs, should be an excellent proxy for 
cluster total mass (Ida Silva et aT1l2lMlMotl et al.ll2005t 
lNagail; 2006: Kr avtsov et al.ll2006h . This presupposes that 
one can determine accurately the value of the integrated 
SZE signal to some mass-scaled cluster radius, as well 
as perform an accu rate calibra t ion o f this relationship 
observationally. In iMotl et all (|2005f ) , we showed that 
the value of Y500, the integrated Compton y parame- 
ter inside a radius where Sp/p — 500 (with respect to 
the critical density) accurately measures the cluster to- 
tal mass inside that same radius. In order to accurately 
calibrate this relationship, one must measure the Y-M 
relationship to high precision at low redshift, where clus- 
ters have combined SZE/X-ray observations to use. Then 
the measured relation, scaled for redshift, can be used to 
determine masses of SZE-selected clusters, which should 
be indentifiable to high redshifts, in order to constrain 
cosmology. 

Since SZE radial profiles are of relatively poor qual- 
ity so far, a direct determination of the compatibility of 
the X-ray /3-model parameters with the SZE profiles in 
clusters can not currently be done observationally. Here 
we compare the values of /3-model parameters for a large 
number of simulated clusters when fitting the X-ray sur- 
face brightness profile, the SZE radial Compton y profile, 
and jointly fitting both X-ray and SZE profiles to a com- 
mon /3-model. We show that the use of X-ray parameter 
values leads to biased estimates of I500, as well as the 
integrated gas mass M$oo tgas . 



We discuss our numerical simulations in Section 2, re- 
sults of the analysis in Section 3, the consequences of the 
use of X-ray /3-model parameters for estimating the Y500- 
M500 relationship for clusters in Section 4, and discussion 
and conclusions in Section 5. 

2. NUMERICAL SIMULATIONS 

Our simulations use the hybrid Eule rian adaptive 
mesh refinement hydro/N-b ody code E nzo (|Q'Shea et al.l 
(|2005h ; |http:/ /cosmos.ucsd.edu/enzo [) to evolve both the 
dark matter and baryonic fluid in the clusters, utiliz- 
ing the piecewise parabolic method (PPM) for the hy- 
drodynamics. With up to seven levels of refinement in 
high density regions, we attain spatial resolution up to 
~ 16 ft. -1 kpc in the clusters. We assume a concordance 
ACDM cosmological model with the following parame- 
ters: il h = 0.026, n m = 0.3, n A = 0.7, h = 0.7, and 
(is = 0.9. Refinement of high density regions is per- 
formed as described in IMotl et~aT1 (|2004h . 

We have constructed a catalog of AMR refined clus- 
ters identified in t he simulation volume as described in 
lLoken et al.l ((2002). The catalog of clusters used in this 
study includes the effects of radiative cooling, models 
the loss of low entropy gas to stars, adds a moderate 
amount of supernova feedback due to Type II super- 
novae in the zones where stars form, and is identified 
as the SFF (S t ar F ormation with Feedback) catalog in 
iHallman et al.l (|2006l ). The catalog includes clusters with 
total mass (baryons + dark matter) greater than 10 14 M Q 
out to z=2 in the simulation. This catalog includes 
roughly 100 such clusters at z=0, and has 20 redshift 
intervals of output, corresponding to a total of roughly 
1500 clusters in all redshift bins combined. 

For parts of the analys is, we have cleaned the clus- 
ter sample as described in [Hallman et al.l (|2006f L remov- 
ing cool core clusters and obviously disturbed clusters. 
This is done primarily to more closely mimic obser- 
vational studies of cluster properties. We examine by 
eye all cluster projections and remove any with obvi- 
ous double peaks in the X-ray or SZE surface bright- 
ness images, have disturbed morphology within R = 
XhT 1 Mpc, exhibit edges consistent with shocked gas, 
or those that have cool cores. The cool core clusters are 
identified as those with a projected emission-weighted 
temperature profile which declines at small radius, or 
those with strongly peaked X-ray emission. The cool 
co re clusters are elimin ated because as we have shown 
in IHallman et al.l (|2006h . they lead to strong biases and 
increased scatter in estimates of cluster physical proper- 
ties when standard observational methods are used. We 
show in our previous work that eliminating obviously dis- 
turbed clusters has a minor, but measurable effect on 
observationally derived quantities. Indeed the simple as- 
sumptions that are typically made in the observational 
derivation (e.g., spherical symmetry, hydrostatic equilib- 
rium) clearly do not apply to such clusters. We also use 
two orthogonal projected images for each cluster. The fi- 
nal cleaned sample contains 493 cluster projected images 
from a series of evolutionary epochs from z=0 to z=2. 

3. RESULTS 

3.1. Isothermal (3-model Fits 

We have generated images of our simulated cluster cat- 
alog by projecting the physical quantities on the grid to 
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Fig. 1. — Upper panels: Left panel shows comparison of fitted /3 values for clusters in our numerical sample at z=0 out to r2500 using 
the X-ray and SZE images. Right panel shows comparison of fitted value of r c for same clusters. Middle panels: Same as above, except 
values are for profiles fitted out to r500. Lower panels: Same as above, but for profiles fitted out to T200- 



get the value of the SZE Compton y parameter and the 
X-ray surface brightness. For the X-ray, we have used a 
simple bremsstrahlung emissivity. It is a simpler calcu- 
lation than using a model X-ray emissivity, and we will 
show later that the values of the /3-model parameters are 
nearly identical irrespective of which emissivity calcula- 
tion we use. 

From these images, we have created radial profiles in 
annular bins, which we subsequently fit to /3-model pro- 
files. Each cluster, for both the SZE and X-ray profile, 
has a set of /3-model parameters that describe it, though 
there is some degeneracy in the parameters in each case. 
The three /3-model parameters, namely the normaliza- 
tion (Sxo or j/o ), the value of the core radius (r core ), and 



the power law index (3 are left as free parameters in the 
fit. We have fit to each of three limiting outer radii in 
both the SZE and X-ray case, ^2500, ^500 > ar >d r 200- In 
each case, the subscript indicates the average overdensity 
with respect to critical inside that radius. This radius is 
calculated from the simulation data using the overden- 
sity of the dark matter. We have used the full sample 
for this part of the analysis, but have excluded the cool 
core region (typically «100kpc) of the profiles from the 
fitting procedure. 

We have found that there are significant differences 
in the values of the fitted model parameters depending 
on whether the X-ray or SZE profiles are used. Figure [T] 
shows the comparison of the values of (3 and r core plotted 
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Fig. 2. — Upper: Ratio of SZE /3 fitted to radial profiles from 
projections of full SFF simulated sample of clusters to /3 fitted 
to X-ray radial profiles for the same clusters. Median values and 
la error bars for the ratio fiszE I Px—ray in z=0 clusters fitted 
to each of three radii T2500i r 500i an d ^200 using the isothermal 
/3-model. Lower: Median values and lcr error bars for the ratio 
r C ore,SZE/r C ore, X-ray in z=0 clusters for the three radii r2500, 
r 50 o, and r 2 oo- 

against one another for individual z=0 clusters in our 
catalog for the three limiting outer radii used. It is clear 
that there is a large amount of scatter in the relationship 
for the P values, and also at rsoo and r2oo, a definite 
discrepancy. Fitting SZE profiles results in a consistently 
higher value of (3 than does fitting X-ray profiles. While 
there is scatter in the compared values of r core , there 
is general agreement within the errors. However, since 
there is a degeneracy between r core and j3 in the fitting, 
some scatter is expected. Since the values of (3 and r core 
are used directly in the equation for the density profile, 
inconsistent values for these parameters lead to different 
deduced density profiles, and discrepant values for the 
cluster gas mass. 

Next, looking at Figure^ we show the lcr scatter in the 
ratio of the SZE deduced parameters and the X-ray pa- 
rameters for each of our three fiducial radii for all clusters 
in our catalog. While there is not a statistically signif- 
icant bias in r core , one clearly exists in (3 at the larger 
radii. To answer the question of why the isothermal (3- 
model fit for X-ray and SZE observations give a different 
set of model parameters in each case, consider the simple 
argument in the next section. 

3.2. The Problem with the Isothermal Model 



TABLE 1 

Mean parameters for Universal Temperature Profile 



Parameter 


2 = 


2 = 0.1 


2 = 0.25 


2 = 0.5 


2 = 1.0 


T 


1.25±0.06 a 


1.27±0.06 


1.31±0.08 


1.37±0.09 


1.37±0.14 


S 


0.51±0.21 


0.51±0.19 


0.42±0.13 


0.43±0.12 


0.53±0.22 


a 


1.17±0.46 


1.12±0.40 


0.91±0.34 


0.80±0.28 


0.90±0.45 



The X-ray and SZE surface brightness functions result- 
ing from the isothermal /3-model density profile integra- 
tion have distributions at r > r c which depend on radius 

as 

S x (r) ocr 1 ^ (6) 



and 



y(r) oc r 



1-3/3 



(7) 



We have fit the temperature profile of our simulated 
clusters to a universal temperature profile (UTP) of the 
form 

T{r) = <T) 500 T ( 1 + ( —Y) , (8) 

where T500 indicates the average spectral (in this study, 
emission- weighted) temperature inside r 500 . T , a, and S 
are dimensionless fitted parameters to the spherically av- 
eraged (from the three-dimensional simulated data) tem- 
perature profiles of all clusters at each redshift in the 
simulations used in this study. The mean values of these 
parameters for clusters at select redshifts from our simu- 
lations are shown in Table [TJ We have used the redshift- 
specific mean value for each cluster in the analysis. 

If the cluster's true gas temperature declines with ra- 
dius with the above described dependence 

-26 



T(r) oc r 



(9) 



then the cluster observable profiles have dependence at 

r > r c 

.1-6/3-4 



and 



SjtWocr 1 "^ (10) 

y(r) oc r 1 -^' 26 . (11) 

Our simulations show a typical value of 5 = 0.5 (see Table 
[T|). If we assume the true cluster density profile is a [3- 
model modified in this way by a UTP for the temperature 
profile (and <5=0.5), and that the true value for (3 set 
by the cluster density profile is /3=0.8, then the radial 
dependence of X-ray surface brightness and SZE surface 
brightness are 

S x (r) oc r" 4 - 3 (12) 

and 



y(r) oc r 



-2A 



(13) 



Finally, if we then fit an isothermal /3-model to cluster 
profiles with the above dependence, setting powers equal 
for the X-ray 

1 - 6/3 = -4.3 (14) 



we would get = 0.; 
fit would give us 



and for the SZE, the isothermal 



1-3/3= -2.4 



(15) 



or (3 = 1.13. 

While the declining temperature profile has a relatively 
small effect on the (3 value extracted from an isothermal 
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TABLE 2 

Median Parameter Values from Profile Fitting to T2oo 



surface brightness 



Method 


P 


+1(T 


-la 


Tcore (kpc) 


+ 1(7 


-la 


Iso X-ray 


0.84 


1.02 


0.70 


168 


330 


99 


Iso SZE 


1.05 


1.27 


0.88 


196 


340 


130 


U-/3 X-ray 


0.81 


0.97 


0.67 


165 


320 


96 


U-/3 SZE 


0.82 


0.97 


0.69 


160 


260 


103 


Density 


0.70 


0.87 


0.55 


137 


242 


73 


Iso Raymond-Smith 


0.82 


1.03 


0.67 


168 


318 


96 



fit to the X-ray surface brightness (+10%), there is a 
larger effect on the value of (3 in the SZE case (+41%). 
This is not surprising given the difference in temperature 
dependence between Equations 2 and 4. Table [2] shows 
the median values of fitted (3 in our cluster sample out to 
r 2 oo, with la scatter. Indeed the variation between the 
X-ray and SZE fits is consistent with the simple analy- 
sis shown above. Also, one line of the table shows the 
result of fitting the true density profiles extracted from 
the simulated data. Though the differences between the 
true value and the SZE and X-ray fitted values are big- 
ger than expected from a simple analysis, the true value 
is indeed smaller than that fitted from the emission pro- 
files. There are additional sources of bias which may 
be introduced due to dumpiness in the ICM and devia- 
tions from spherical symmetry which m ay contribute to 
the larger differen ce in fitted parameters (Sulkaiicnl ll999t 
iNagai et aLll2000f) . There is also the de generacy in r c and 
(3 in any model fitting to these profiles which can con- 
tribute to differences. We also show in the last line of the 
table the result of fitting the isothermal /3-model to the 
profile of X-ray surface brightne ss resulting from a pro- 
jection of the Raymond-Smith (Brickho use et al1l!995f ) 
model emissivity. It is clear that the statistical proper- 
ties of the fitted parameters are nearly indistinguishable 
from those fitted to the simple bremsstrahlung model. 

It is therefore clear that the effect of a declining tem- 
perature profile in real clusters is to alter the fitted values 
of (3 in this simple model, such that the SZE isothermal 
(3 is larger than both the X-ray, and the "true" value 
associated with the density profile of the cluster. Indeed 
a similar analysis bv lAmeglio et all (|2006l ) has shown re- 
sults consistent with ours in this regard. We also find 
no significant trend in the SZE/X-ray parameter ratios 
with redshift. 



3.3. Fitting to a UTP Modified Model 

Since the variation in model parameters between X- 
ray and SZE fitting appears to be due to the failure to 
account for the radial dependence of temperature in the 
intracluster medium, it makes sense to use a model that 
includes this radial dependence for fitting the surface 
brightness. Our non-isothermal /3-model for the surface 
brightness is created by integrating the expression for the 
X-ray surface brightness and Compton y parameter us- 
ing the standard /3-model for the density, and the UTP 
for the temperature. We refer to this model as the U-/3. 
Integrating Equations 2 and 4 with the substitution of 
the /3-model density profile (Eq. 1) and the UTP (Eq. 
8) results in the fitting relations for the X-ray and SZE 




respectively, where Ix{b) and Isz(b) are line integrals 
described in the Appendix. 

As described in Section 3.2, the parameters a and S are 
fixed in the fitting to the surface brightness distributions, 
and result from fitting of the spherically averaged (from 
the three-dimensional simulated data) temperature pro- 
files from the simulated clusters. The average values of 
the parameters in each redshift bin are used for fitting 
the surface brightness profiles of a cluster in that same 
redshift bin. When fitting the cluster profiles to these 
relations, we find that there is now consistency in the 
values of (3 and r c between X-ray and SZE fitting, shown 
in Figure [3] While there is still some scatter, the offset 
is removed. Table [2] shows the results for all clusters in 
the sample when fitting to an isothermal /3-model or a 
non-isothermal (the U-/3) model. The value of (3 is vir- 
tually identical in SZE and X-ray fitting when using the 
U-/3. Additionally, they result in similar distributions of 
r c values. 

Note that this model has the same number of free pa- 
rameters (three) used to fit the surface brightness as a 
standard /3-model, the others are fixed from simulations. 
We should note the well-known caveat about /3-model 
fitting, that the values of r core and (3 are somewhat de- 
generate in this fitting. 

4. CONSEQUENCES OF THE (3 
INCOMPATIBILITY: EXAMPLE CALCULATION 
OF Ysoo VS M 

Using the isothermal /3-model parameters from a fit 
dominated by the X-ray data to do SZE or combined 
X-ray/SZE analysis introduces an additional error or 
bias to various derived quantities. Here, we character- 
ize the nature of these errors on the measured Y500 vs 



M: 



500, gas 



relation for clusters. As described in the intro- 



duction, it has been shown that the integrated Compton 
y par ameter inside rsoo i s an excellent proxy for total 
mass (Ida Silva et al.l f2000l iMotl et al.ll2005l : iNagall [20061 : 
iKravtsov et al.ll2006l ). We examine whether one can ac- 
curately determine the true value of I500 (and M^oo,gas) 
from a cluster observation using standard observational 
techniques. 

Generating a value for Y500 using a /3-model fit dom- 
inated by the X-ray emission, one must extrapolate to 
rsoo using the model for the SZE emission. Currently 
existing cluster SZE profiles from interferometric instru- 
ments do not constrain individual /3-model parameters 
well. Since we have shown that the isothermal /3-model 
paramters for SZE and X-ray cluster profiles are differ- 
ent, there is an error introduced. Additionally, there are 
errors introduced in the calculation of the cluster gas 
mass, Msoo^gas, since in any individual cluster, the (3- 
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Fig. 3. — Upper: Ratio of SZE f3 fitted to radial profiles from 
projections of full SFF simulated sample of clusters to j3 fitted to 
X-ray radial profiles for the same clusters. Median values and lc 
error bars for the ratio 0s ze/ Px—ray m z= clusters when fitted 
to each of three radii r2500, t"500j and f200 when fitting with a 
U-/3 model. Lower: Median values and lc error bars for the ratio 
r C ore,SZE/r C ore, X-ray in z=0 clusters for the three radii r 2 5oo, 
r 50 o, and r 2 oo- 

model for the density deduced from the observations is 
only approximately correct. 




Fig. 4. — Upper: Integrated Compton y parameter inside rsoo 
from simulated clusters at z=0 epoch plotted against integrated 
mass inside the same radius. Stars are for the true values on 
the simulation grid, red diamonds represent values deduced by fit- 
ting to isothermal /3-models, blue triangles are for fitting to a U-/3 
model. Lower: Log deviation of isothermal and U-/3 points from 
the best-fit scaling relation for the true values. 

by extrapolating the model for y(r) using the parame- 
ters P and r c from the X-ray fitting, and the value of yo 
generated by fitting the SZE profile with (3 and r c fixed. 
M500, gas is estimated by integrating the f3 model for den- 
sity out to an overdensity of 500. This method effectively 
measures t^qq from the fitted gas density profile in order 
to get the integrated gas mass. For the U-/3 cases, the val- 
ues of a and S are fixed at the mean values for the whole 
sample of clusters at each redshift, recognizing that these 
values would be provided by simulations, and would not 
be left as free parameters in the observational analysis. 



4.1. Estimation of Y and M 

We have used the simulated cluster sample described 
above (SFF) to determine the systematic errors intro- 
duced in this procedure. We have analyzed the simulated 
clusters such that our synthetic observations are analo- 
gous to unbiased, high signal-to- noise observations of real 
clusters. The simulated observations are idealized, since 
no instrumental effects or foreground/background source 
removal are simulated. We have performed our calcula- 
tions on projected SZE and X-ray images generated from 
each simulated cluster, assuming the X-ray profiles could 
be determined to rsoo- 

Using this method, and assuming that the values of (3 
and r c are defined by fitting of the X-ray surface bright- 
ness profile, we have characterized errors by taking the 
model parameters from the X-ray fitting and calculated 
the estimated values of Y500 and M5oo,gas from each clus- 
ter. We then compare the value to the true value taken 
from the simulation grid. 

We have used the isothermal model fits and also the 
fits to the U-/3 models. In each case, Y500 is determined 



4.2. Comparison of Isothermal and U-/3 Methods 

On first inspection of FigureSJ it is not clear that either 
the isothermal or U-/3 generates a more accurate result. 
The plot shows values from all clusters from the z=0 
simulation epoch. The true values of F500 and Msoo, ff as 
are plotted as stars, the isothermal model values as dia- 
monds, and the U-/3 model values as triangles. 

The plot of deviation of these values from the best fit 
to the true values shows some difference between these 
two methods. However, it is important to remember that 
there are deviations from both the true integrated SZE, 
and the true gas mass, so these points deviate in both di- 
mensions of the plot. Since errors in both are correlated, 
due to the dependence on the determination of rsoo from 
the /3-model, this should reduce the apparent separation 
of points from the best fit, even when the values of both 
parameters are in error. 

This important point is illustrated in Figure El which 
shows the the ratios of estimated to true values of Y500 
and Msoo^gas for each cluster. In both the isothermal 
(upper plot) and U-/3 (lower plot) methods, there is a 
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Fig. 5. — Upper: Plot of the ratio of estimated Y500 to its 
true value (y-axis) versus the ratio of estimated Msoo.gas to its 
true value for 493 simulated cluster projections for isothermal 
/3-model fits. Lower: Ratio of estimated Y500 to its true value 
versus the ratio of estimated Msoo,gas to its true value for 
U-/3-model fits. Lines indicate location of perfectly accurate 
estimation of each quantity. 



clear trend in over and underestimation in both values. 
These figures also illustrate nicely the improvement in 
values one gets with the \J-(3 method. The U-/3 points are 
clustered more strongly around the correct values than 
are the isothermal estimates. This effect is illustrated in 
Figure [H where we have plotted the true values of y and 
M, with lines indicating the corresponding estimates of 
those values for the isothermal and XJ-/3 methods. There 
are typically overestimates in both estimated Y500 and 
-^500, gas, and the use of the U-/3 typically brings the 
estimated values back toward the true. 

The outliers in these plots are also interesting. The 
clusters which give strong overestimates of Y and M gas 
typically result from line-of-sight overlap of multiple 
structures. This enhances the integrated SZE signal in 
projection, but overestimates the true value for the main 
cluster. The clusters with low estimated values of Y 
and/or M gas typically have poor quality /3-model fits (as 
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Fig. 6. — Plot of eight clusters on the Y500--M500, gas relation 
chosen to span a wide mass range, but otherwise randomly. Stars 
indicate the true values from the simulation, red diamonds the 
isothermal /3-model estimates, and blue triangles the U-/3 model 
estimates. Lines link points associated with the same cluster. 

measured by a x 2 statistic) , and appear to have system- 
atically low fitted values for the core radius. Some of 
these clusters appear to have weak cool cores, not meet- 
ing the cool core criteria, and thus were not excluded 
from the sample. It is interesting that the clusters with 
mass underestimates do not all generate underestimates 
of Y. This effect warrants further study. 

Table [3] also shows the result of this analysis. When 
using isothermal models, the median estimated value of 
Y500 is 23% larger than the true value (the mean is a 
28% overestimate). Additionally, when estimating the 
gas mass, the median est imated value is a 16% overesti- 
mate. In previous work (Hall man et alJ [2006). we have 
shown that the magnitude of the overestimate via X-ray 
/3-model methods is lower than this, but in that case we 
assumed one could correctly calculate the value of rsoo 
with no error. In either case, this overestimate is con- 
sistent with our p r evious work, and with that of others 
(|Mohr et alj Il999t iMathiesen et all I1999D . The overes- 
timate in mass results from substructure, merging and 
other physical processes not described by a simple model. 
In the current analysis, the value of rsoo is deduced di- 
rectly from the /3-model. In contrast, the U-/3 method 
gets closer to the correct value for Y500, with a median 
value 11% higher than the true value. For the mass, the 
U-/3 method gets an estimate close to the true mass, a 
median overestimate of only 4%. 

Additionally, using the U-/3 method reduces the scatter 
in Y500 values, as shown in Figure [7] The distribution 
of values is more sharply peaked, with a smaller high 
end tail in addition to a reduced bias compared to the 
isothermal method. In mass, the result is a reduced bias 
in the median values with a slight improvement in the 
scatter as shown in Figure [5J 

5. CONCLUSIONS 

There is an inconsistency between X-ray and SZE fitted 
model parameters that leads to a bias in deduced values 
of Y§oo and M$oo t g as when using isothermal /3-models. 



TABLE 3 

Ratio of Estimated to True Values of Simulated Cluster 
Properties to rsoo 





Isothermal 


v-p 


(2/500, est/2/500, true) 


1.28 


1.13 


(2/500, est/2/500, true), median 


1.23 


1.11 


la upper 


1.50 


1.24 


la lower 


1.08 


1.03 


{^500, gas, est /M500, gas, true) 


1.17 


1.05 


(-W500, gas, est /M50O, gas, true), median 


1.16 


1.04 


la upper 


1.31 


1.18 


la lower 


1.06 


0.94 
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Fig. 7. — Histogram of 493 simulated cluster projections of ratio 
of estimated I500 to true value for isothermal (black lines) and U-/3 
(red lines) methods. Dotted lines indicate median values for each 
distribution. 



Using our U-/3 model reduces the bias and scatter, re- 
sulting in a more precise and accurate determination of 
the y-M scaling relation for clusters. 

X-ray and SZE radial profiles from the same galaxy 
cluster can not be fit well by identical isothermal /?- 
models. The stronger dependence of the SZE emission 
on temperature leads to problems when the cluster tem- 
perature declines with radius. In contrast, the strong 
dependence of the X-ray emission on density minimizes 
the error introduced by variations in cluster temperature. 
We show that fitting either the X-ray or SZE profiles to a 
modified, non-isothermal /3-model corrected by the inclu- 
sion of a univeral temperature profile for clusters results 
in better consistency with the real values of Y500 and 

-^500, gas- 

We expect this inconsistency should affect measure- 
ments of the Hubble constant at some leve l, indeed it ap- 
pears to result in a bias in Hq , as shown bv lAmeglio et ail 
(2006). There will certainly be an effect on gas fraction 
determinations, though it will also depend on any bias in- 
herent in X-ray hydrostatic total mass estimates. While 
current precision in SZE/X-ray derived quantities is not 
high enough to reveal effects at the 10-20% level of the 
type described here, we expect it will be higher with 
larger samples and new instruments available in the near 
term. When contemplating clusters as precision cosmo- 
logical tools, effects as this level must be considered. 

While some ha ve used double /3-models 
(|LaRoque et al.l 120061 ) to fit cluster radial profiles 
with some success, a model like the U-/3 is preferred, 
since it is physically motivated by the observed ICM 
properties. This work illustrates a simple modification of 
the isothermal /3-model to more realistically account for 
cluster physics, and thus more accurately and precisely 
measure cluster properties. 
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Fig. 8. — Histogram of 493 simulated cluster projections of ratio 
of estimated Afsoo, gas to true value for isothermal (black lines) and 
V-/3 (red lines) methods. Dotted lines indicate median values for 
each distribution. 
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APPENDIX 

DERIVATION OF UTP SURFACE BRIGHTNESS MODEL 

To derive a surface brightness model resulting from the UTP modified /3-model, we take the standard /3-model for 
density, 

/ , x 2\ 

n e {r)=n e0 1+ - , (Al) 



r c 



and the UTP model for temperature, 



V ar 500 

and substitute them into the integral of the X-ray surface brightness, 



T(r) = <r)Boor { 1+f—J I , (A2) 



Sx(b) = 47r(1 X +z)4 J n e (r)n H (r)A(T(r))dl. (A3) 
In the simple bremsstrahlung case, A(T) = AoT 1 / 2 . In that case, 

&<»>-&.(l +(£)") (;£;)') <A4) 

where 

2n e o 2 Ao((r) 500 ro) 1 / 2 

5X0 = 4™(1 + ,) 4 ' (A5) 

where k = (j,h/ fJ>e, the ratio of the mean molecular weights of hydrogen and electrons. For r c < arsoo, 

r c ( ( b \ 2 \ 1/2 „ {I „„ S 1\ „ /{ 1 „„ 5 r c 2 + 6 2 



where i<2,i is Gauss' hypergeometric function and -B is the Beta function defined by 

B(x,y) = IffiO ;a!j y >0. (A7) 
r(a: + y) 

For r c > ar 500 , 

<*« - ^ (' + (^) 2 ) " B {*» + H F - (* b 3 « 4 1 - ' « A8 » 

Similarly for the SZE, where 

V = f cp^n (r)-^d^, (A9) 
J m e cr 

the substitution results in 

3/3 _r 

2/(6)=yo ( 1+ (^) 2 ) ( 1+ (q^) 2 ) Isz{b)l (A10) 

and 

y = 2o- T n e o 5 • (All) 

m e cr 



For the two cases of r c < arsoo an d r c > arsoo respectively. 

^(l+(^) 2 )^B(K^+S-^F it (8^..,. 2 
2 I \rj j ^2 2 2) 2A \ '2' 2 ' a 2 r 500 



Wi) = 7 ( 1+ ( -) I g(o,-f + g-o) J , a,i(^oS-^ + ^ 1 - -2 2^1.2 )■ ( A12 ) 



and 

r„„™ - fi j. , , L „ , .. . 

\ar 500 J J V 2 ' 2 2 J V 2 ' 2 ' 2 ' r c 2 + 6 2 

Figure [9] shows the values for each of the line integrals as a function of projected radius normalized to rsoo for typical 
model parameters. We have created a calcul ator for the line in tegral in the Interactive Data Language (IDL), which 



is available on-line at http://solo. Colorado. edu~hallman/UTP/Ib_calc. tar. gz 
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